Strain-dependent Splitting of Double Resonance Raman Scattering Band in Graphene 
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Under homogeneous uniaxial strains, the Raman 2D band of graphene involving two-phonon 
double-resonance scattering processes splits into two peaks and they altogether redshift strongly 
depending on the direction and magnitude of the strain. Through polarized micro- Raman measure- 
ments and first-principles calculations, the effects are shown to originate from significant changes in 
resonant conditions owing to both the distorted Dirac cones and anisotropic modifications of phonon 
dispersion under uniaxial strains. Quantitative agreements between the calculation and experiment 
enable us to determine the dominant double- resonance Raman scattering path, thereby answering 
a fundamental question concerning this key experimental analyzing tool for graphitic systems. 
PACS numbers: 63.22.Rc, 63.20.dk, 73.22.Pr. 



The effects of external mechanical perturbations on 
physical properties of graphene are attracting much at- 
tention because of the possible realization of synthetic 
electromagnetic fields [lHfl and determinations of its fun- 
damental material parameters Since electronic 
structures and phonon dispersion are modified under 
applied strains [lH3j, 5 16], the Raman spectrum - an 
important diagnostic tool for graphitic systems - will 
show significant variations. Recent experiments demon- 
strate that the Raman G band redshifts and splits into 
two peaks under strain because of symmetry breaking 
of the doubly degenerate Ei g phonons [7H9(. Further- 
more, one can determine the Griineisen parameter of 
graphene 0JM and identify its crystallographic ori- 
entation [7|, [8fl. 

Unlike the Raman G band, the scattering process of 
the Raman 2D band involves electronic states and TO 
phonons near the K and K' points of the Brillouin 
zone so that the strain-induced anisotropy of the elec- 
tronic band structure @, Q and the phonon dispersion 
111 |l5[ must manifest themselves in the Raman spectra 
of strained graphene. Although there have been experi- 
mental analyses assuming isotropic TO phonon softening 
(l6| and independent theoretical studies regarding the 2D 
band of strained graphene 15| , a comprehensive and sys- 
tematic study considering changes in both electronic and 
phonon structures is still lacking. Moreover, because sev- 
eral resonant scattering processes contribute to the 2D 
band, a fundamental question concerning the dominant 
double- resonance process remains to be resolved (l7| . 

In this Letter, we present a comprehensive analysis 
of the changes in electronic energy bands and phonon 
dispersion of a single-layer graphene under homogeneous 
uniaxial strains by combining polarized Raman measure- 
ments with an analysis based on first-principles calcula- 
tions and determine the dominant scattering path of the 
double-resonance Raman scattering process. As the mag- 
nitude of the strain increases, the Raman 2D band is split 
into two peaks, both of which redshift. Moreover, two 



distinct strains applied along armchair and zigzag crys- 
tallographic directions are identified and the frequency 
shift rate for each split Raman 2D peak is strongly depen- 
dent on the strain direction. From theoretical analysis, 
we demonstrate that the anisotropic TO phonon soften- 
ing together with distortions of Dirac cones is a dominant 
factor responsible for the observed effects. Furthermore, 
the polarization dependence of the relative intensities of 
the split 2D band components reveals contributions of 
different resonant scattering paths, thereby establishing 
fundamental understanding of the double-resonance Ra- 
man scattering process in graphene. 

Single-layer graphene samples were prepared on acrylic 
substrates with 50 x 10 x 1.3 mm 3 dimensions by using the 
micromechanical cleaving method from natural graphite 
flakes [lU . Single-layer graphene samples were identified 
with micro- Raman spectroscopy 1JJ-21|. The results 
from two samples with special orientations will be com- 
pared. Strain was applied by bending the substrate with 
a specially designed jig. The Raman spectra were ob- 
tained using a polarized micro-Raman system using the 
514.5-nm line of an Ar ion laser as the excitation. Other 
experimental details have been previously published [22| . 

As the magnitude of uniaxial strain (e) increases, the 
G band redshifts and splits into twopeaks, G~ and G + 
(not shown) , as was reported earlier pHSI . The shift rates 
are du G -/de = -33.4 cm-7% (-33.0 cm-7%) and 
duj G+ /de = -14.5 cm-7% (-12.9 cm-7%) for sample 
A (B). They are essentially the same for the two sam- 
ples as expected for small strain from symmetry and are 
in agreement with Ref. 7 but larger than the values in 
Ref. 8, probably due to difference in strain calibration. 
We also obtain the Griineisen parameter of 2.2 ± 0.1 and 
the shear deformation potential of 0.93±0.04, in excellent 
agreement with previous estimations 0, [23| . The angle 
between the strain direction and the zigzag direction of 
the graphene lattice, ips, was found to be 34.9 ± 0.2° 
(52.7 ± 0.5°) for sample A (B) from the polarization de- 
pendence of the relative intensities of the G~ and G + 
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peaks @, HI- Since cps = 30° and 60° for strain applied 
exactly along the armchair and zigzag directions, we will 
henceforth refer to samples A and B as A-strain and Z- 
strain samples, respectively. 

When strain is applied, the 2D band splits into two 
peaks which redshift as the strain increases [Figs. HJa) 
and (b)]. Unlike the G band, the frequency shift rates of 
the 2D~ and 2D + peaks of the A- and Z-strain sam- 
ples are significantly different from each other: for A 
strain, dio 2D -/de = —63.1 cm~ 1 /% and duj 2 o+/de = 
—44.1 cm _1 /%, whereas for Z strain, duj 2D -jdt = 
-67.8 cm-7% and duj 2D+ /de = -26.0 cm^/Yo. Again, 
these values are larger than those in Ref. 16, presumably 
due to difference in strain calibration. 

The dependence of strain-induced 2D band splitting 
on the strain direction offers a unique opportunity to 
examine the strain-induced anisotropy of the electronic 
and phonon bands. The 2D band comes from the four- 
step Stokes-Stokes double-resonance Raman scattering 
as illustrated in Fig. EJa) [H, HE 0]. Theoretical 



calculations 26, 1271) suggested that the scattering pro- 
cesses involving the smallest momentum transfer (inner 
process) and the largest momentum transfer (outer pro- 
cess) are dominant contributions, but it is still not clear 
which of the two is the dominant one [13] • For a given 
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FIG. 1. (Color online) Evolutions of the 2D bands of (a) 
A-strain and (b) Z-strain samples as a function of uniaxial 
strain. In (a), the 2D~ and 2D + peaks are obtained with 
9- m = 6 OVL i — 0° and f? m = OVL i = 90°, respectively. In 
(b), the spectra are measured with f? m = 8 OVL ^ = 50° and 
deconvoluted into two Lorentzian peaks. The positions of 
the 2D~ and 2D + peaks of (c) the A-strain and (d) Z-strain 
samples as a function of strain. The solid lines are linear 
fits to the data. The inset shows the polarization geometry, 
where f? m , on t, and ips are the angles that the incident laser 
polarization, the analyzer axis, and the zigzag direction make 
with respect to the strain axis, respectively. 



laser wavelength, the momentum of the emitted phonon 
is determined by the electronic band structure and the 
phonon dispersion near the K and K' points. For un- 
strained graphene, the scattering processes involving the 
three K' points around a given K point [denoted by 1, 2 
and 3 in Fig. HJb)] are completely equivalent, and hence 
the 2D band appears as a single peak. When A strain is 
applied, the reciprocal lattice is distorted, as in Fig.^c); 
one of the three K' points moves away whereas the other 
two K' points move closer to the K point. Therefore, 
the two types of scattering processes involve phonons 
with different momenta, resulting in a splitting of the 
2D band. For Z strain, the distortion of the reciprocal 
lattice is reversed as in Fig. Eld). In Figs. EJc) and (d), 
the scattering process 1 involves phonons with momenta 
in the DKS direction, whereas processes 2 and 3 involve 
phonons with momenta in the FRM direction. 

In order to analyze the observed splitting of the 2D 
band quantitatively, modifications of both the electronic 
band structure and the phonon dispersion due to strain 
must be taken into account. Strain shifts Dirac points 
away from the K or R points and tilts and distorts the 
Dirac cone so that the group velocity depends on the di- 
rection in the Brillouin zone 0, ■ Strain also modifies 
the phonon dispersion. The observed splitting and soft- 
ening of the 2D band thus result from a convolution of 
the electronic band structure and the phonon dispersion 
modifications. So, we performed calculations on elec- 
tronic band structures of strained single-layer graphene 




FIG. 2. (Color) (a) Double-resonance Raman scattering pro- 
cess. Inner (q m ) and outer (g ou t) processes are indicated. 
Reciprocal lattice diagram for (b) unstrained, (c) A-strained, 
and (d) Z-strained graphene, showing strain-induced distor- 
tions. The K' points are now designated as R points, and 
there are inequivalent high symmetry points M and S, the 
midpoint between the K and R points. 
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FIG. 3. (Color) Calculated phonon momenta involved in the 
(a) inner (q m ) and (b) outer (<? ou t) processes for the scatter- 
ing in the FKS and FRM directions and A and Z strains. 
Phonon dispersions near the K or R points for 2% strain ap- 
plied in the (c) armchair and (d) zigzag direction. Filled and 
open stars are for the 2D band in unstrained graphene, filled 
and open squares for the 2D + peak, and filled and open cir- 
cles for the 2D~ peak. Calculated strain dependences of the 
2D peaks for (e) inner and (f) outer processes. 



based on the first-principles self-consistent pseudopoten- 
tial method Q using the generalized gradient approx- 
imation for exchange-correlation functional p9| and on 
their phonon dispersions by using density-functional per- 
turbation theory 28|, 30] . The ion core of carbon atoms is 
described by an ultrasoft pseudopotential [3l| . A fc-point 
sampling of 48 x 48 x 1 grid uniformly distributed in the 
two-dimensional Brillouin zone is used in self-consistent 
calculations and a 6 x 6 x 1 grid is used to calculate the 
dynamical matrices. The obtained electronic structures 
and phonon dispersions for graphene without strain are 
in good agreement with other studies [optical phonon 
frequencies at r and K points are 1581 cm" 
1295 cm" 1 



(E 2g ) and 

{A[), respectively] [32 1. 

The calculations show that the dominant contribution 
to the observed splitting and strain-direction-dependent 



frequency shifts originates from anisotropic changes of 
TO phonon branches with distorted Dirac cones as pre- 
sented below. Since modifications of the electronic struc- 
ture will change the resonant conditions as discussed 
above, the scattered phonon momentum will change their 
magnitude depending on their directions that are de- 
termined by momentum conservations under the strain. 
Figures [3t[a) and (b) are the calculated phonon momenta 
satisfying the resonant conditions for the inner (<7j n ) and 
outer (<7 ut) processes, respectively. They are signifi- 
cantly different for the rKS and rRM directions and 
for A and Z strains. Figures [3jc) and (d) illustrate the 
modified phonon dispersions near the K or R points for 
2% strain applied in the armchair and zigzag directions, 
respectively. In the figure, the corresponding resonant 
frequencies for the Raman 2D band are indicated by filled 
(inner) and open (outer) symbols, respectively. It should 
be noted that the phonon dispersions along the rKS 
and rRM directions are significantly different, especially 
away from the K point. This is in direct contradiction 
to the assumption used in Ref. 16 that the phonon soft- 
ening rate is orientation independent. Figures Efe) and 
(f) summarize the strain dependences of the 2D peaks 
for the inner and outer processes, respectively. It is clear 
that the outer process is not consistent with our exper- 
imental data shown in Fig. [TJ A linear fit to the calcu- 
lated values up to 2% of strain gives, for the inner process 
[Fig. He)], duj 2D -/de = -70cm-y% and du 2D+ /de = 
— 43cm~ 1 /% for A strain and du> 2 o-/de — — 66cm~ 1 /% 
and doj 2 £>+/de = — 24cm~ 1 /% for Z strain, in excel- 
lent quantitative agreement with the experimental data. 
Therefore, our Raman data for strained graphene clearly 
demonstrate that the inner process is the dominant one 
in the double- resonance Raman scattering. We also 
note that the softening of the TO phonon is more or less 
isotropic at high symmetric points (K and R) as assumed 
in Ref. 16, but the strain phonon momentum satisfying 
double-resonance conditions deviates from those points 
where the effects of anisotropic softening arc significant. 
It should be noted that the 2D splitting rates for A and 
Z strain could not be explained in Ref. 16 with a model 
that assumed an isotropic TO phonon softening. 

Finally, the intensities of the 2D~ and 2D + peaks de- 
pend strongly on the polarization direction of the inci- 
dent laser [Figs. QJa) and (b)], corroborating our anal- 
ysis. The intensities of the 2D~ and 2D + peaks were 
measured as a function of the incident polarization an- 
gle (# m ) in steps of 10°. The analyzer is kept parallel 
to the incident polarization, 6- m = OVl ^, which prefer- 
entially selects phonons in the direction orthogonal to 
& n 22, [33|. In A-strained graphene, for example, the 



phonons involved in the scattering processes of 2 and 3 in 
Fig. [2^c) have the same frequency, whereas the phonons 
for the process of 1 differ. Since each of the three scatter- 
ing processes contributes to the 2D band, the peak corre- 
sponding to the process of 2 and 3 should have a contri- 
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FIG. 4. (Color online) Evolutions of the 2D band of the 
(a) A-strain and (b) Z-strain samples as a function of the 
incident laser polarization relative to the strain axis (# m ), 
under 0.97% strain. The analyzer is parallel to the incident 
laser polarization. Polar plots of the 2D~ and 2D + bands of 
the (c) A-strain and (d) Z-strain samples as a function of # m . 



bution twice that of 1. This is demonstrated in Figs.^Jc) 
and (d). In Fig. SJc), the 2D~ band has one sinusoidal 



component fitted well to I 2 d- °c cos 



61), whereas 



the 2D + band is fitted to two sinusoidal components of 



2n/3) + cos- 
2 = 3.1° [Hi 



9 in - 2 - 47T/3), 
In Fig. Iljd), the 



I 2D + oc cos 4 (6» m - (j> 2 
where (f>i = 5.1° and 
2D~ band has two sinusoidal components fitted well to 
I 2 d- °c cos 4 (f9 m — 4>\ - 7tt/6) + cos 4 (# m — 4>i — ll7r/6), 
whereas the 2D + band, with one sinusoidal component, 
is fitted to I 2 d+ °c cos 4 ((5 m — 2 — tt/2), where 4>i — —6.0° 
and 4> 2 = -7.5°. 

In conclusion, the strain-induced splitting and redshift 
of the Raman 2D band are found to depend on the direc- 
tion of the applied strain with respect to crystallographic 
orientation. Comparison of experimental data with first- 
principles calculations shows that anisotropic modifica- 
tions of the phonon dispersion together with changes in 
electronic structures are their origins. Furthermore, the 
dominant inner scattering process is demonstrated to re- 
solve a controversy regarding the nature of the Raman 
2D band. 

Note added: Recently, we became aware of the publica- 
tion of related work on the Raman 2D band of strained 
graphene 16[, which was reported simultaneously with 
ours recently 35!]. 
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